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УДК 539.32 
А.О. ВАТУЛЬЯН, А.Ю. БУРЬЯН, А.В. ОСИПОВ 


ОБ ИДЕНТИФИКАЦИИ ПЕРЕМЕННОЙ ЖЕСТКОСТИ 
ПРИ АНАЛИЗЕ ПОПЕРЕЧНЫХ КОЛЕБАНИЙ БАЛКИ 


Представлен метод нахождения амплитудной характеристики для консольно закрепленной балки перемен- 
ной жесткости на основе анализа интегрального уравнения Фредгольма 2-го рода, предложен метод опреде- 
ления закона изменения жесткости на основе решения задачи Коши. 

Ключевые слова: балка переменной жесткости, коэффициентная обратная задача, интегральное уравне- 
ние Фредгольма 2-го рода, сплайн-аппроксимация. 


Введение. Коэффициентные обратные задачи теории упругости — относительно новый и значи- 
мый раздел математической теории упругости, применяемый в теории неоднородных покрытий, 
фундаментостроении, сейсмологии, био- и наномеханике и других областях. Для этого класса за- 
дач можно выделить модели, в которых идентифицируемые дифференциальные операторы име- 
ют постоянные коэффициенты (линейная теория анизотропного упругого тела, линейная теория 
вязкоупругого тела при дифференциальной форме определяющих уравнений), и модели, в кото- 
рых требуется идентифицировать неоднородные свойства (например, в геофизике при описании 
деформирования литосферных плит и при разведке полезных ископаемых, в горной механике при 
анализе напряженного состояния в окрестности выработок и особенно предварительного напря- 
женного состояния при изучении наноразмерных объектов, в биомеханике при исследовании раз- 
личных тканей и вибрационных воздействий на них в целях идентификации) [1]. 

Методы определения жесткости важны для процедуры идентификации объектов в различ- 
ных областях естествознания. Главная проблема при исследовании задач подобного типа - это 
формулировка операторной связи между искомыми коэффициентами дифференциальных опера- 
торов и известными (измеренными) функциональными зависимостями. Представлены результаты 
исследования обратной задачи для дифференциального оператора 4-го порядка, описывающего 
колебания балки переменной жесткости. 

Постановка задачи. Рассмотрим установившиеся колебания упругой балки длины / переменной 
жесткости Е] с частотой ‹. Уравнение имеет вид [2]: 


(Ех, )м”(х, ))" — реЕихх, ) = 0, (1) 
где Е - модуль Юнга; / — момент инерции; р -— плотность; Е’ - площадь поперечного сечения. 
Считаем, что балка на конце х, =0 жестко закреплена, а на конце х, =/ - действует на- 


грузка. Рассмотрим два варианта нагружения на конце — 
задача 1: нагружение сосредоточенным моментом; 
задача 2: нагружение сосредоточенной силой. 
Соответствующие граничные условия имеют вид — 








задача 1: 
и(0) = и” (0) =0; ЕЛ -м"(=М ; (ЕЦх,)-м"(х „1 =0; 
задача 2: 
м0) = и’ (0) =0; ЕК." =0; (Ех) и" (х,))' мы Р. (2) 
Приведем краевую задачу (1), (2) к безразмерному виду, введя безразмерную координату 
2 
х=2^1, а также безразмерные параметры и функции: В(х) = Е: ‚К = а 
/ ЕК ри 
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Таким образом, краевая задача (1), (2) в безразмерных переменных имеет следующий 
вид: 
(Вхуи’(х))" - Кикх)=0. (3) 
Поскольку решение задач 1, 2 пропорционально нагрузке, то примем следующие модель- 
ные граничные условия — 
задача 1: —и(0)=и’(0)=0; В(1)-м’(Т) =1; 


задача2: —и(0)=и’(0)=0; В(1):м’(=0; 


Основная цель исследования: 
1) разработка метода построения функции и(х,К) при различных значениях параметра 





[А 








(4) 


К при произвольных законах изменения жесткости В(х) (монотонных, немонотонных, кусочно- 
непрерывных); 

2) разработка метода построения решения обратной задачи, которая состоит в 
нахождении функции В(х) по информации о функции и(х,^) при некотором значении д, ; 

3) проведение серии вычислительных экспериментов, позволяющих оценить 
эффективность предлагаемых методов. 
Исследование прямой задачи. Прямая задача состоит в нахождении функции "›, удовлетво- 
ряющей краевым условиям при заданном законе изменения В(х) в некотором диапазоне измене- 


ния А. Поскольку коэффициенты дифференциальных операторов переменны, то в аналитическом 
виде найти решение нельзя, возможно только численное решение задачи. Чаще для этого ис- 
пользуют метод конечных разностей или метод сведения к интегральному уравнению Фредгольма 
2-го рода [3]. 

Предлагая способ сведения краевой задачи к интегральному уравнению Фредгольма 2-го 
рода, опишем процедуру сведения на примере задачи 2. 

Для построения интегрального уравнения введем функцию: 

и(х) = В(х)м"(х). 

Далее перейдем от краевой задачи (3), (4) к системе дифференциальных уравнений 2-го 

порядка относительно функций и(х) и и(х): 


и"(х)- К‘ иХх) = 0, 


5 и(х 5 
О-о (5) 
В(х) 
Соответствующие граничные условия приобретают вид 
их0) = и”(0) =0, 
(6) 


и =0, и’) =1. 
Интегрируя уравнения (5) дважды по переменной х, находим представления функций 
и(х) и и(х) в следующей форме: 


и(х) = ес —5и(5)46 + Ах + 4, 
(7) 


их) = | (х-5) Г ЧЕ+В,+В,. 


Постоянные 4,, 4., В,, В, находим из граничных условий (6): 
1 1 1 1 
В =8,=0, А =1- А" [9(4Б, 4, = [94 -1-#' [а-Вжба к" [5 94Е-1. 
0 0 0 0 
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Таким образом, имеем соотношения для определения искомых функций: 


их) = [«-9 г Е, 





их =" [(<- Эх х[оав+и" [баЕ-1=(х-0+# [*-Э\94&- 


= [с -904Е = 2-9-# |«-Вн(4& 
0 х 
Объединив эти представления, получим: 


д- [9 т и — | (5). (8) 


Поменяем порядок интегрирования в двойном интеграле в представлении (8), переписав 
его в виде интегрального уравнения Фредгольма 2-го рода: 





1 
их) = А* | К(х, 5) и (5) + Г, (9) 
0 
причем его ядро непрерывно для широкого класса функций В(х) и представимо в форме 


етве-В аа 


кое 89 
[-ве-ю Е ав 
В© 
10-9&-0 
и 


Для задачи 1 соответствующее уравнение получается аналогично и имеет вид (9), но с 
другой правой частью, а именно: 


2-9 Ч. 
В(5) 

Уравнения вида (9) можно исследовать либо методом последовательных приближений, 
либо численно. Наиболее распространенным методом является метод коллокаций. Дискретиза- 
цией на равномерной сетке й =/М№, х, =Ай, А=0,...,№, хх =1, заменой интеграла квадратур- 
ной формулой Симпсона [5] и использованием метода коллокаций уравнение (9) приводится к 
исследуемой численно системе линейных алгебраических уравнений: 


го) = Г. 


М№ 
и(х,)- УТ А,К(х,х,)и(х,) = 1(х,), =Ъ2,..,М. (10) 


Анализ конечномерного оператора (10) позволяет эффективно строить численное реше- 
ние и определять как резонансные частоты, так и амплитудно-частотные характеристики для ба- 
лок переменной жесткости. 

Для оценки точности предложенной схемы вычислений сравним численное решение, по- 
лученное в результате использования составленной программы, и точное решение, полученное 
аналитически при В =1, которое выражается через функции Крылова [2]. 
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Расчеты показывают, что относительная погрешность между этими решениями не пре- 
восходит 0.5% в диапазоне частот до первого резонанса, что свидетельствует о достаточно эф- 


фективной вычислительной схеме. 
Результаты расчетов амплитудно-частотной характеристики представлены в виде 


графиков. 
На рис.1 результаты решения задачи 2 изображены сплошной линией при А =5, а ромби- 


ками — при А =1 для В(х)=1+х”. 


%(х) 


—0,1 


-0,2 


—9,3 





Рис.1. Решение прямой задачи 2 


На рис.2 изображен результат решения задачи 1, сплошной линией показано решение при 
К =5, а ромбиками — при К =1 для В(х)=1+ 2 (10лх). 


\(х) 
0,4 
0,3 


0,2 


0,1 


-0,1 


-0,2 





Рис.2. Решение прямой задачи 1 
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Обратная задача. Целью решения обратной задачи является восстановление функции В(х) по 
известным в наборе точек значениям функции и_(х,А‚). При этом рассмотрены варианты поста- 
новки обратной задачи: 

задача 3: Восстановление функции В(х) по значениям функции и(х,А,) (которая в рам- 
ках вычислительного эксперимента взята из решения прямой задачи 1); 

задача 4: Восстановление функции В(х) по значениям функции "Хх,А,) (которая в рам- 
ках вычислительного эксперимента взята из решения прямой задачи 2). 

Заметим, что подобные обратные задачи рассмотрены и в другой постановке [6, 7]. В 
данном случае дифференциальное уравнение (3) относительно функции В(х) имеет 2-й порядок 
и для ее определения необходимо решать задачу Коши. Некорректность обратной задачи в этой 
постановке состоит в том, что необходимо находить вторые производные от функции, заданной в 
наборе точек. Поэтому на лервом этапе аппроксимируем функцию э_х) сплайнами 5-й степени 
[5], что позволяет находить вторую производную с достаточной степенью точности. На втором 
этапе построим операторное соотношение, связывающее функции В(х) и иХх). Осуществим это 
построение на примере задачи 4. Дважды интегрируя уравнение (3) от | до х, получаем: 


х 


[Водь"соу а = [ (водо) -1= &* [5(6\Е, 


и далее имеем: 


х-1-* [©«-Б и 


В(х)= 11 
(х) и "О (11) 
Для задачи 3 аналогичное соотношение запишется в виде 
1 
1-4" [(х- ЭЙ 
В(х) = и 12 
(х) и") (12) 
Дискретизация позволяет найти значения искомой функции: 
1 
1 
Ви) = м1 [с - ЭМ |. (13) 
и (%) ы. 


1 


Для вычисления определенного интеграла [2 — Е) (5 а 6 заменим подынтеграль- 
ХК 
ную функцию на каждом отрезке разбиения многочленом Лагранжа 1-й степени с узлами 1 и х 
[5]. Это соответствует замене на каждом отрезке разбиения кривой на секущую. В результате 
получаем известную формулу трапеций: 


[© - В) иХЕМЕ = п, = Жен )+ (Хр — Хи Ух )+ „Е 


№ 


+ (ХЕ — Ху Ух )+ В Ух И 
Применение этой формулы к (11) дает расчетные зависимости. 


829 


Физико-математические науки 








В(*,)= м НТ наиб )+ (к хо ина) 


1+... (№ Руда) 69 — хн их ]] : 


Путем аналогичных построений получаем расчетные формулы для задачи 3: 
1 4 
В(х = " я | > [а (© жи )+ 
и (+) 
(= Жи хьь2 + ЗН хуи + 12 (% = ху) х» )) 
Как видно из представления (11), задача 4 хуже обусловленная, нежели задача 3, так как 
на конце х=1 имеет место неопределенность, поэтому в этой точке необходимо либо задавать 


значение функции В(х), либо экстраполировать функцию из области х<!1. Этому факту есть 


простое объяснение в рамках рассмотрения задачи Коши. Для задачи 3 граничные условия (4) 
дают два необходимых условия в задаче Коши, а для задачи 4 — только одно, поскольку из усло- 
вия В()-и’(Т) =0 и положительности В(х) следует м”(1Т) = 0, и второе условие пропадает. 


На основе описанного метода проведены эксперименты по реконструкции различных ви- 
дов неоднородностей: гладкой монотонно возрастающей функции, гладкой монотонно убываю- 
щей функции, гладкой немонотонной функции, функции, имеющей на отрезке [0, 1] несколько 
стационарных точек, кусочно-разрывной функции. 

На графиках представлены вычислительные эксперименты для различных видов неоднородности. 


Рис.3 иллюстрирует результат задачи 4 по восстановлению функции В(х)=1+х? (гладкой 
монотонно возрастающей функции) при А =1: 

м(х) 

2 


1,5 


0,5 


: 





0 0,2 0,4 0,6 0,8 1 


Рис.3. Решение обратной задачи 4 


В качестве примера на рис.4 изображен результат восстановления функции 
В(х)= 1+ зт? (10лх) (функции, имеющей несколько стационарных точек) в задаче 3 с входными 


данными, зашумленными аддитивным образом с амплитудой 8 =10^. 
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У(х) 
2 у >> > > © 
1,5 
1 ®> $ ® ®* ® ® ® > х 
0,5 
0 0,2 0,4 0,6 0,8 1х 


Рис.4. Решение обратной задачи 3 с зашумленными входными данными 


В случае разрывных законов неоднородности проведено разработанный метод усовершен- 
ствован: входной массив данных исследуется на предмет точек разрыва (рассматривается изме- 
нение первой и второй производной функции Хх)). При обнаружении точек разрыва в их окре- 
стности проводится экстраполяция функции, что улучшает качество восстановления. 

На рис.5 изображен результат восстановления функции: 


1 при безе не. 
В(х) = : - 4 к =0,9. 
5 при —<х<-, 
3 4 


5 Фохоофофоь® > 





0 0,2 0,4 0,6 0,8 1х 
Рис. 5. Решение обратной задачи 4 для кусочно-разрывной функции 
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Выводы. Проведена серия вычислительных экспериментов по реконструкции переменной жест- 
кости различных типов. Представлена методика решения коэффициентной обратной задачи при 
анализе поперечных колебаний балки переменной жесткости. Разработан эффективный метод 
решения при точных и зашумленных входных данных. 
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А.О. УАТУТУАМ, А.У. ВУКУАМ, А.\. ОЗТРОУ 


ОМ ТРЕМТТЕТСАТТОМ ОРГ УАКТАВЕЕ ВТСТОТТУ ОМОЕВ 
ТНЕ АМАЕУ$1$ ОЕ ТНЕ ВЕАМ ТАТЕКАЕ ОЗСТЕЕАТТОМ 


Тре теод оЁ атрЖиае сПагацей$@с езбтаНоп Гог {Пе сапШеуег Беап оР {Ле уацаЫе пааКу оп {Ле 
Баз!5 ог {Ле апа\$!5 ог РгедПот и\едга! едцаНоп оР {Ле зесопа Кта 15 ргеземе4. ТАе деегттавоп 
теёоа ог поаКу сПапдта оп {Не Баз оЁ Саиспу зедиепсе апау$5. 

Кеу мгога$: Беат оЁ {Пе уайаЫе поаКу; соеаепЕ тшуегзе ргоМет; РгедПот и\{еодга! едиаНоп оЁ 
{Пе есопа Кта; эрИтпе арргохитаНоп. 
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